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^ ' Abstract 



We examine the scaling properties of an evolving network of strings in Minkowski 
spacetime and study the evolution of length scales in terms of a 3-scale model pro- 
posed by Austin, Copeland and Kibble (ACK). We find good qualitative and some 
quantitative agreement between the model and our simulations. We also investigate 
small-scale structure by altering the minimum allowed size for loop production E c . 
Certain quantities depend significantly on this parameter: for example the scal- 
ing density can vary by a factor of two or more with increasing E c . Small-scale 
structure as defined by ACK disappears if no restrictions are placed on loop pro- 
duction, and the fractal dimension of the string changes smoothly from 2 to 1 as the 
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resolution scale is decreased. Loops are nearly all produced at the lattice cut-off. 
We suggest that the lattice cut-off should be interpreted as corresponding to the 
string width, and that in a real network loops are actually produced with this size. 
This leads to a radically different string scenario, with particle production rather 
than gravitational radiation being the dominant mode of energy dissipation. At 
the very least, a better understanding of the discretisation effects in all simulations 
of cosmic strings is called for. 

1 Introduction 

Cosmic strings formed during a cosmological phase transition at an energy scale of about 
10 16 GeV are sufficiently massive to seed structure formation. However, calculations of 
the gravitational effects from strings are hampered by uncertainty over the statistics of 
the evolving network [[l], 

From lattice simulations of string formation ||, a picture has developed of a string 
network at formation consisting of a scale-invariant distribution of closed loops together 
with a percentage of infinite brownian strings crossing the Universe. The subsequent 
evolution would consist of the growth of the step length of the brownian infinite string 
network as the strings attempt to straighten out, and energy loss via loop production 
during reconnection. 

Early analytic work identified the key property of scaling where at least the gross 
properties of the network can be characterised by a length scale, roughly the persistence 
length or the interstring distance £, which grows with the horizon ||]. This result was 
supported by subsequent numerical work ||. However, further investigation revealed 
dynamical processes, including loop production, at scales much smaller than £ || 0]. In 
response, Austin et al developed a model describing the network in terms of three scales 
|| (hereafter ACK). These scales are: the usual energy density scale £, a correlation 
length £ along the string and a scale ( relating to local structure on the string. It seemed 
likely from the ACK model that £ and £ would scale with ( growing slowly, if at all, until 
gravitational radiation effects became important when £/£ « 10~ 4 ||. Progress has also 
been made using alternative analytical approaches in |I| |Tl|] . 

Understanding how the network scales is important because it simplifies the process 
of building a model representing the evolving network. Predictions of the power in the 
anisotropies in the Microwave Background require the two-time correlation functions for 
various components of the string energy- momentum tensor evolving over O(10 3 — 10 4 ) 
expansion times. As present simulations cannot achieve this range, constructing a model 
is necessary. We presented a simple model of the two-time correlation functions in a 
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non-expanding background in |T2 



In this paper we present the results of a numerical study of evolving strings in a 
Minkowski space-time. We demonstrate the scaling behaviour of correlation functions 
of p and q vectors (dynamical quantities which are linear combinations of tangent and 
velocity vectors) along the string. These functions may be used to complement analytical 
studies || and also in existing formalism to predict small-scale anisotropy in the cosmic 
microwave background [13]|. The largest length scale associated with these correlations, 
£, and the interstring distance £ are both significantly smaller than the causal horizon 
which suggests that the effects of the expanding background on the correlation functions 
are small. Following Coulson et al [14|, we think of the Minkowski network of strings 
as modelling more realistic strings, by mapping Minkowski time to conformal time and 
Minkowski space to comoving coordinates. 

As a test for the ACK model, we examine the rate equations for the evolving length 
scales developed in ACK as they apply to a Minkowski space-time with no gravitational 
radiation, and compare the results with our simulations. Such a comparison is problem- 
atic because there are a number of parameters which are difficult to calculate directly. 
In particular, the parameter k which controls the behaviour of ( is hard to measure. We 
have showed good quantitative agreement between the ACK model and our simulations 
for the approach to scaling of £, and qualitative agreement for (. 

A major problem is understanding the relationship between small-scale structure 
along the string and loop production. Bennett and Bouchet || provide evidence that 
strings possess a fractal substructure arising from kink production during reconnection. 
This "intermediate" fractal, constant at fixed time, spreads over a larger range of scales 
as the network evolves and by the end of the simulation covers a range between the 
resolution scale and the persistence length £. Most of the loop production occurs at 
scales near the resolution scale. Shellard and Allen || also observe loop production near 
the scale of resolution and an "intermediate" fractal in the long string, although they 
attribute the fractal to the effect of initial conditions and a restrictive dynamical range. 

The build up of small-scale structure is allowed by the Austin et al analysis: their 
original guesses for the model parameters yield the increase in small-scale structure ob- 
served in the expanding Universe codes. However, it is not well known that their model 
can also predict that the small-scale structure (as measured by () is absent. 

Indeed, we find that if there are no restrictions on loop production, then small-scale 
structure will disappear: all three length scales defined in ACK will have the same 
magnitude. We also find that the region of constant fractal dimension, which we observe 
if we artificially constrain loop production, disappears. This has interesting implications 
for the energy loss mechanism from cosmic strings. 
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The standard scenario is that any build up of small-scale structure will only be checked 
by the back reaction from the string's own gravitational field, which only becomes effective 
when £/£ is small (of order 10~ 4 ) [|§. Loop production will then grow with the horizon, 
albeit at a scale much smaller than £. Although much smaller than £, these loops are 
vastly bigger than the string width and will be topologically stable, losing energy through 
gravitational radiation. 

If there is no small-scale (as measured by () and loop production occurs at the smallest 
possible scale - presumably the string width - then energy loss will be dominated by 
production of GUT quanta fT5(| . Furthermore, particle production from cuspy regions 
may occur without any conventional loop production at all. This will obviously allow 
the string scenario to avoid any gravity wave bounds QTJ] , but the implications for the 
decay products of the string quanta are less clear. Recent calculations of the flux of high 
energy decay products from string quanta have been compared with observed fluxes to 
put an upper bound on the proportion of string energy going into particles [|16|]. In the 
standard scenario these bounds can be avoided fl7| . However, if energy loss is dominated 
by particle production such bounds will constrain the GUT physics if the GUT scale 
string scenario is to survive. 



2 Algorithm and Code 

In Minkowski space the string equations of motion are 

X" = X, (1) 

where X = dX/drj and X' = <9X/<9s; rj is time and s is a space-like parameter along 
the string. The spatially flat FRW metric is conformal to the Minkowski metric and 
we regard the Minkowski time as conformal time, in the limit that the expansion of the 
Universe goes to zero. 

X = X(s,?7) is a position three-vector which satisfies the constraints 

X' • X = 0, (2) 

X' 2 + X 2 = 1. (3) 

The first constraint ensures no tangential velocities, and @ ensures that the energy of a 
segment of string is proportional to its length when measured in units of s. 



Using a development of a code written by one of us previously [0, 18[, we solve the 



wave Eq. ([!]) with the Smith- Vilenkin algorithm [I9||. This algorithm uses the exact 
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finite difference solution to ([!]), 

X(s, r] + 5) = X(s + 5, rf) + X(s - 5, r,) - X(s, rj-S). (4) 

If the string points are initially defined on the sites of a cubic lattice (N5) 3 , then 
(|) ensures that they remain on the lattice at time steps of 5. The discretised gauge 
conditions require that for these initial conditions, string links are restricted to three 
types: stationary (|X| = 0, |X'| = 1), diagonal (|X| = l/y/2, |X'| = l/\/2) and cusps 
(|X| = 1, |X'| = 0). 

Because the string points lie on the sites of the lattice, identifying crossing events is 
simple. When two strings cross, they intercommute with a probability which is set to 
Pi. For most of this work we set Pi = 1, subject to the condition that reconnection 
does not create a loop smaller than the threshold E c . Loops with energy greater than 
or equal to a threshold value of E c are allowed to leave the network, while reconnections 
are forbidden for loops with energy equal to E c . Forbidding reconnections allows energy 
to leave the network fairly efficiently; otherwise it takes much longer for the effect of the 
initial conditions to wear off. This feature may also model more realistic networks as, in 
an FRW Universe, small loops will decouple from the expansion and are highly unlikely 
to reconnect. 

An interesting value for E c is the minimum segment length 25. The loops produced 
in this case are "cusps": confined to one lattice site and travelling at the speed of light. 
The cutoff then becomes a lattice cutoff. 

The ability to alter E c allows us a certain amount of control over small-scale structure 
and related quantities. For example, as E c is increased from 25 to 85, the rms velocity 
of a segment along the string increases from 0.36 (very close to the matter era value in 
Bennett and Bouchet) to 0.46 (close to the radiation era value). This corresponds to a 
build up of fast-moving small-scale structure. 

Initial string configurations are generated using the Vachaspati-Vilenkin algorithm 0, 
which mimics the breaking of a U(l) symmetry during a cosmological phase transition. 
In this approach, phases from the minimally discretised U(l) group are placed randomly 
on the sites of a cubic lattice. Strings (or anti-strings) are identified as passing through 
a face with a net winding about the manifold. The strings are then joined up within the 
lattice cells to from a network (any ambiguity if two strings leave and two strings enter a 
lattice cell is settled randomly). The domains, and consequently the initial network step 
size £oj are °f constant size. This initial configuration is defined on a cubic lattice with 
fundamental lattice spacing 5. We can alter the initial correlation length £o in terms of 5. 
If each initial segment length is made up of a large number of individual stationary string 
links, peculiarities arise in the scaling solution as entire segments can be annihilated into 
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Figure 1: Strings formed with the Vachaspati-Vilenkin algorithm in a (128(f) 3 box with 
£o = 85. This snapshot is taken shortly after formation. 



loops of energy E c moving at the speed of light. Consequently we add structure in the 
form of cusps to break up long straight segments. Cusps are string links confined to one 
lattice site which move at the speed of light. Figures [l] and § show a Vachaspati-Vilenkin 
string network shortly after formation and after a period of evolution. We have studied 
simulations on lattices ranging from (64<5) 3 to (10245) 3 . We varied £o from 85 to 645 and 
E c from 25 to 85. Periodic boundary conditions are imposed throughout. We restrict the 
evolution time to half the box size, as after this time causal influences have propagated 
around the box. 

An ensemble typically consists of 30-50 runs. Unless otherwise stated, measurements 
were taken from networks created in a (1285) 3 box with £o = 85. 

3 Left and right movers 

The most general solution to Eq. ([I]) for cosmic strings in a Minkowski spacetime is 



This solution may be considered as made up of "left-moving" (a) and "right-moving" (b) 
pieces. We will consider the dynamical quatities p and q defined 



(5) 



p = b' = X + X', 



(6) 
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Figure 2: The strings pictured in the previous figure after the evolution period with 
E c = A5. 




Figure 3: The strings pictured in Figure [I] after the evolution period with E c = 25. The 
number of small loops is suppressed. 
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Figure 4: Scaling C pp (upper curve) and C pq correlation functions for £ = 86 and E c = 26. 

q = -a' = X - X', (7) 

where p 2 = q 2 = 1. The solution (|5|) may then be represented by a pair of curves for 
p and q on a unit ( "Kibble- Turok" ) sphere. String intercommutation can excise a loop 
between s\ and s 2 if 

P ds(p - q) = 0. (8) 

At such an intercommutation event, kinks are created as abrupt changes in the p and 
q vectors. Left and right moving kinks travel away from the intercommutation site 
in opposite directions. Correlations between ps and qs along the string are extremely 



important [§, [L3|] as they relate to loop production, but cannot easily be calculated a 
priori. If s = Si — s 2 

C PP (s,r]') = (p(si) ■ p(s 2 ))\ v=v >, (9) 

C pq (s,7]') = (p(si) ■ q(s 2 )}\ v=v >, (10) 

These two correlation functions contain all the information needed as C qq (s) = C pp (s) 
and C qp (s) = C pq (—s). We find that these functions approximately scale with £ i.e. they 
are functions of s/£ only. We calculated these correlation functions over an ensemble 
of 50 realisations for various E c . The results for the scaling functions for E c = 26 are 
plotted in Figure @. The ensemble averages are fitted well by the simple functions: 

C PP (\s\/0 = (1 - uo)e- a ^ + ue- b ^, (11) 

C pq (s/0=de-^/^ 2 . (12) 



The first of these was used in |J for demonstration purposes, but is in fact a good 
approximation to our measurements. It is particularly good for the E c = 25 case. The 
gaussian form for C pq becomes less good away from the peak as the tail may become 
exponential. The gaussian fit is however within the ensemble errors. The parameters for 



parameter 


E c = 25 


E c = 45 


E c = 85 




0.57 ± 0.05 


0.45 ± 0.03 


0.25 ± 0.03 


a 


5.4 ± 1.2 


20 ± 3 




b 


0.54 ± 0.1 


0.35 ± 0.04 


0.17 ± 0.03 


d 


-0.34 ± 0.02 


-0.20 ± 0.02 


-0.08± 0.02 


f 


0.55 ± 0.1 


0.30 ± 0.05 


0.03 ± 0.01 


e 


0.70 ± 0.1 


1.1 ± 0.1 


3.8 ± 0.2 



Table 1: Parameters for model C nn and C nn functions. 



the model functions are given in Table |I[ The parameters a and b pick out two length 
scales along the string. The effect of increasing E c , and therefore small scale structure, 
is to decrease the small scale 1/a and increase the large scale 1/6, relative to £. However, 
for E c > 45 the small scale set by 1/a is no longer constant and C pp and C pq do not truly 
scale. 

The related correlations between the tangent and velocity vectors along the string are 
shown in Figure [5]. These functions can be calculated as simple sums of ( |Tl~D and ( |T2"D 
using equations (|) and (0). 



4 Length scales and rate equations 

Analytic work as so far led to three dynamical length scales characterising the network. 
Austin, Copeland and Kibble M have considered network evolution in terms of the dy- 
namics of left- and right-moving kinks and have defined three length scales used in their 
analysis which we can calculate using our simulations: 

(i) £, the familiar energy density length scale 

(ii) £, a persistence length along the string 

(iii) (, a measure of the small-scale kinky structure. Their definitions are 

£(rj) = Jn/pi{v), (13) 
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Figure 5: The scaling functions for correlations between the following vectors: tangent- 
tangent (— ), velocity- velocity ( ) and tangent- velocity (— •) for £ = 85 and E c = 26. 

oo 

m = fds(p( Sl ).p(s 2 )), (14) 



= ((l-u)/a + u/b)£, (15) 

<T\V) = -V^(P(0).P( S ))| S =0, (16) 
= {{l-u)a + ub)(-\ (17) 

where £ and £ are also expressed in terms of the specific correlation function given in 

(O). 

The angular brackets () indicate averaging over the long string, which is defined as 
all string having energy greater than £, and over an ensemble of realisations. 

ACK developed rate equations for the three length scales. From their (expanding 
universe) analysis it seemed most likely that the network would enter a transient scaling 
regime, in which £ and £ scaled but ( did not. (£ would scale only after gravitational back 
reaction was taken into account, when C/f ~ 1CT 4 ) . In order to study this solution, we 
considered these rate equations and made an approximation suitable to our simulations: 
we ignore terms involving Hubble expansion or gravitational radiation. Then the rate 
equations become: 

i - (i8> 

<i9) 
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Figure 6: Plot of £ and £ for a (1285) 3 box with £ = 45 and £ c = 25. 



(20) 



In Eq. (|IBD , c characterises loop production. ACK consider it as function of ratios of 
length scales, although it should tend to a constant in the scaling regime of our simula- 
tions. In Eq. (|19"D, x is a geometrical factor affecting the frequency of intercommutation 
between uncorrelated string segments; u is a parameter in the correlation function C pp 
as described by Eq. (]TT|) . The function / also relates to loop production. In Eq. (|20|). 
k relates to the efficiency of small-scale structure removing itself from the long string 
network. We change variables to x = £/r], x = C,/r] and z — Q/r] to express the above 
rate equations: 

(21) 

1/ J- - J: (22) 

//: ■— —+'^-z. (23) 



cx 

T]X = — - X, 

2x 

-XX 2 I 
ujx 2 2 



x. 



x 



5 The Approach to Scaling 

As in ACK we express each length scale as the product of time and a scaling variable 
so that £ = xrj, £ = xt] and ( = zrj. It has been widely reported || ||, [7| that £ scales. 
Inspection of Figure ^| indicates the possibility that £ and £ are not in fact scaling, but 
growing with a power of rj less than 1. A correction to the standard scaling scenario is also 
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suggested by a recent analytical calculation A study of the exponent n in r] n as the 
network evolves shows that it is indeed tending towards 1 throughout the simulations. 
At late times we find no significant departure from the scaling value n = 1. For a 
(1285) 3 box, n = 0.96 ±0.05 which may indicate than £ is still on an approach to scaling. 
However, as the box size is increased, n becomes ever closer to 1 with n = 1.000 ± 0.003 
for a (10245) 3 box. The errors on £ are larger, but similar conclusions apply. To extract 
the scaling values for the two length scales, we fit the data C/rj towards the end of the 
simulation to x* + where and x* are free. This is motivated in part by the 

approach to scaling predicted from the ACK model at the end of this section. 

The resulting scaling values x* and x* are given in Tables |] and || The figures are 
given with 1-sigma error bars from averaging over an ensemble of realisations. 



Co 


E c = 25 


E c = 45 


E c = 85 


85 


0.174 ± 0.018 


0.130 ± 0.012 


0.089 ± 0.006 


165 


0.170 ± 0.010 


0.141 ± 0.010 


0.109 ± 0.004 


325 


0.171 ± 0.011 


0.147 ± 0.006 


0.117 ± 0.006 


645 


0.173 ± 0.011 


0.144 ± 0.003 


0.122 ±0.014 



Table 2: Scaling values for x = d^/drj. 



Co 


E c = 25 


E C = A5 


E c = 85 


85 


0.19 ± 0.04 


0.13 ± 0.03 


0.070 ± 0.03 


1Q5 


0.20 ± 0.04 


0.17 ± 0.03 


0.10 ± 0.03 


325 


0.19 ± 0.03 


0.17 ± 0.03 


0.12 ± 0.04 


645 


0.25 ± 0.04 


0.21 ± 0.04 


0.14 ±0.03 



Table 3: Scaling values for x = d^,/dr]. 



The scaling parameter x* is independent of £o> within ensemble errors. It is however 
strongly dependent on E c to the extent that, with an increase of E c from 25 to 85, the 
scaling density increases by a factor of 2. The scaling values of x* have bigger ensemble 
errors, and this may hide some dependence on £ - However, the results so far show that, 
within the ensemble errors the only dependence is on E c . 

One of the problems testing the model is the difficulty in calculating the parameters 
directly. We must use the model to calculate some of the parameters and see whether 
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this gives consistent results. Assuming that c, /, 
following expressions are stable fixed points. 



J .1, 



c 
2' 



X and k achieve scaling values then the 



(24) 



\ 2u{I 



(2k - l)x 2 Jx 




if 2k - 1 > 
if 2A; — 1 < 0. 



(25) 
(26) 



We are able to calculate c directly from the scaling value of £*. The value of x can be 
measured by counting the number of long string intercommutation events over each time 
step. Following ACK, the probability of a segment of length £ intercommuting in a time 
step 6 is given by: 



e 



(27) 



We assume that the number of segments in the volume V is given by VY(£ 2 £), thus the 
number of long string intercommutating events in a volume V over 5 is given by: 



X V5 



(28) 



Figure [F] shows the number of long string intercommutations as a function of £ for E c = 25, 
together with a slope giving an £ -4 dependence. From the y-intercept we can estimate 
X at 0.05±0.015 for both E c = 26 and E c = 45. We are now able to make a direct test 
of the rate equations by considering the approach to scaling of £. 

Expanding about the fixed points gives (to first order) a matrix equation equivalent 
toEqs. ©-(H): 



where 




(29) 



(30) 



and 



/— - 1 



— 



M 



5™ 

—JT- 



2x 
-2x55 



-kcz* 






+ 



(31) 



kc 



1/ 
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Figure 7: The number of long string intercommutation events (Nlsi) in the simulation 
box at each timestep plotted against £. 



Using the scaling value for x* and given above M reduces to 








—x, 




X* 






i ^4 




LOX% 




I 




— kcz* 




x'i 










(32) 



The eigenvalues are then — 1, 



-|2fc- 1|/ 

and — 1 2 A: — 1|. Thus we expect in the approach to 



scaling (for large rj) the time dependences: 



x(rj) = C\—r\ 
x 



a 



uxi 



rf + X* 



x{rf) = Cif]^ 1 + C 2 rf + x*, 
2(77) = C 3 r] x + z*, 



where 



A 



-\2k - II. 



(33) 

(34) 
(35) 

(36) 
(37) 



To compare with the simulations, we measured x — ^/rj. As predicted from (|53"D, we 
measured a 77 _1 dependence for large 77, although this arises partly from £o/?7- We subtract 
the x*+ar7 _1 dependence from a plot of x, leaving a term proportional to rj e . The resulting 
log- log plot for a simulation with £ = 45 and _E C = 25 is shown in Figure (§), and it gives 
a good single power law fit to e = —1.4. The ACK model predicts e = —1.5 ± 0.2, from 
Eq. (^) and our measured values for x*, and u. Similarly for E c = 45 we measure 
e = — 1.8 ± 0.2 and the model predicts e = —1.7 ± 0.2. We see no evidence that another 
scale enters into the approach to scaling. 
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-2.8 



-3.2 
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In ii 

Figure 8: Log plot of the if term in Eq. (|33f) . 

6 Small Scale Structure 

Austin e£ al characterise small-scale structure using p and q vectors to define a length 
scale, as in Eq. ([H]), from which we get: 

r 1 (^) = (p(0).p'(0)). (38) 

Figure |9] shows the dependence of ( on E c . An E c of 25 allows structure to rapidly leave 
the network and ( to grow. For E c = 85, structure builds up on the string faster then 
loop production can remove it and ( decays. For E c = 25 and E c = 45, ( appears to 
approach a scaling value. For a box of size (128<5) 3 , the exponent n in ( = + zt] n 
achieves n = 0.98 ± 0.03 towards the end of the simulation and tends towards the 
values given in Table |], given for a range of string densities. We note that for both these 
values of E c , £/£ ^> 10~ 4 : much bigger than the ratio required for gravitational radiation 
to become significant. Indeed, for the smallest value, ( scales at the same magnitude as £. 
As we increase E c however, we find that ( tends towards zero growth. This signifies the 
build up of small-scale structure which would eventually trigger effective gravitational 
back-reaction. We would like to encompass the dependence of £ on E c in terms of the 
three-scale model. The ACK model for (, Eq. (p0[), consists of two processes. The first 
is the effect of long string intercommutation. This will serve to create kinky structure 
and tends to decrease (. The second term involves the effect of small loop production. 
It is controlled by a parameter k which reflects the extent to which loop production can 
remove small scale structure. If we assume a scaling regime for £ and £ we can integrate 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 9: The behaviour of ( for E c = 25 (-), E c = 45 ( ) and E c = 85 (— ). The 

plots correspond to different box sizes to ensure that £ / ' E c > 4. 





E c = 25 


E c = 45 


85 


0.11 ±0.01 


0.007 ±0.002 


1Q5 


0.12 ±0.01 


0.008 ±0.002 


325 


0.12 ±0.01 


0.005 ±0.001 


Q45 


0.12 ±0.01 


0.005 ±0.002 



Table 4: Scaling values for z* = d(/drj. 



Eq. fl23|) to get 

(2k - l)ri 2k - 1 , . 

where C is a constant of integration. The evolution of ( is then crucially dependent 
on k. For k > C nas a positive gradient and z tends towards the scaling value 
2* = {2k — l)xl/x- For fc < |, £ has a negative gradient and tends towards 0. From 
Figure |9] we see that both cases occur in our simulations, depending on the value of E c . 
The three-scale model is certainly capable of explaining our simulations, if we consider 
E c as being related to k. This should not surprise us, as increasing E c artificially keeps 
structure on the string that would otherwise leave the network. Unfortunately, we cannot 
calculate k a priori. We can only assume that the rate equations from the three-scale 
model are broadly correct and extract a value for k numerically, by using the scaling 
expression for z*. This gives k = 0.59 (E c = 25), k = 0.51 (E c = 45) and k < \ 
(E c > 65). 
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Figure 10: How the scaling values and z* vary with the intercommutating probability. 



The effect of changing the intercommutating probability on z* is broadly similar to 
the effect of [0]. We compare the two in Figure [II]. The value for z* changes slowly 
until Pj fa 0.35 when it falls rapidly. After this point ( is probably not scaling. 

Another way to characterise small-scale structure along the string is through the 
relationship 

E = R d , (40) 

where E is the average energy between two points separated by a physical distance R. 
The exponent d is scale dependent. On large scales it is shown in || [7j that d tends 
to 2, as the strings become random walks. On small-scales, the string becomes straight 
on average, and d tends to 1. One might expect a smooth interpolation between these 
two regimes, but Bennet and Bouchet note that an "intermediate" fractal d\ (i.e. 
a constant or very slowly varying d) exists on scales between the numerical cutoff and 
below £. They explain this as the build up of small scale structure on over the scales of 
constant fractal. 

Allen and Shellard || observe a similar structure, although they attribute it to arti- 
ficial initial conditions and lack of dynamical range. We observe that this intermediate 
fractal exists for energy cut-offs of 85 and above whereas for E c = 25 it has disappeared 
completely leaving a smooth transition from d = 1 to d = 2. This is shown in Figure |TTJ, 
which is plotted for comparison with similar figures in |J. However, it is not clear that 
small scale structure always reveals itself through a constant d s . As a comparison, we 
calculated E/R from the tangent-tangent correlation function (X'(z) • X'(0)) with the 
parameters in Table [1] for E c = 25. For these model parameters, we find that £ w 0.4£: 



certainly not a small scale. We plot the result in Figure |12] along with the measured 
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Figure 11: Log plot showing the fractal structure of the long string for energy cutoffs of 
E c = 25 (— ), E c = 45 (--) and E c = 85 (-). 

quantity. The agreement is good up to scales well in excess of £. (The discrepancy at 
larger scales may be due in part to model function underestimating the pq correlations 
at large £, as the gaussian form of C pq breaks down.) However, for small scales the 
agreement is good and this indicates a genuine lack of small-scale structure. 

7 Conclusions 

We have demonstrated the scaling behaviour of correlation functions of p and q vectors, 
which are essential in analytic studies of string networks. 

The parameters of model functions for the pp and qq correlation functions C pp and 
C pq , depend on how much small-scale structure in present on the string, which is related 
to the value of the minimum allowed loop energy E c . As E c is increased, small-scale 
structure is allowed to build up and the scaling density can increase by a factor of two 
more. 

We analysed our simulations in terms of the three-scale model proposed by Austin et 
al, suitably simplified for Minkowski space with no gravitational radiation. 

The model for £ describes well the approach to scaling and predicts the two exponents 
for this approach that we measure in our simulations. The model for the "small" scale £ 
is qualitatively consistent with our simulations with critical value k determining whether 
or not ( will scale. 

For the smallest possible value of E c , small-scale structure disappears, as all the scales 
£, £ and £ of the 3-scale model of ACK have the same order of magnitude. Furthermore, 
from the fractal analysis along the string for E c = 25, the transient fractal region dis- 
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Figure 12: Log plot showing the fractal structure with a cutoff of E c = 25 ( — ) together 
with the same quantity calculated from the model tangent-tangent correlation function. 



appears, leaving a smooth transition from straight strings on small-scales to the random 
walk on large scales. Thus the intermediate fractal region seen in the expanding Universe 
codes seems to be an effect of the limit on the size of loops. 

In all cases, loops are predominantly produced at the smallest possible scale E c , 
even when it is set, by the underlying lattice, to 25. Thus, in common with all other 
simulations (apart from those of Albrecht and Turok ||), there is no evidence for scaling 
in loop production. It is therefore time to ask the question: is this the true physical 
situation? It is certainly possible that loops are produced with the smallest possible 
physical scale, which is the string width, and we see from our simulations that every 
other measure of the network has a single scale of order £. Thus there is no real conflict 
with the string scaling hypothesis. 

One objection to the use of the Smith- Vilenkin algorithm is its restrictive set of 
possible values of p and q for the initial conditions (although the subsequent evolution 
is exact). As discussed by Albrecht pi] , this may exaggerate small loop production. 
However, expanding Universe codes have no such restrictions and there is no published 
evidence that loop production happens at any scale other than the loop size cut-off. 
Furthermore, we have found no evidence that the lattice scale enters the scaling solution. 
If there is spurious cusp production (for E c = 25) through a chance configuration of 
p and q vectors such that Eq. (|8|) is true, then one would not expect a stable scaling 
solution. Albrecht suggests that there may also be spurious loop production through 



"back-tracking" |2l| which occurs when individual discretised links or groups of links 



point in a direction opposite to that of the whole segment of length £. In this case, one 
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would expect the lattice scale 5 to affect the evolution: for example, the rate of energy 
loss from back-tracking would most likely increase with £/8 as the chances for a given 
link or group of links to back-track in a time £ increases. 

We are confident that our simulations do approach (and for large simulations reach) a 
stable scaling solution, and from Table § we see that the value of the scaling parameter x* 
is independent of the range of £ over which the simulation is run. We infer from this that 
the spurious processes mentioned above do not critically enter the observed dynamics. 

The prospect that the loops are produced with such tiny sizes is a radical one. It 
means that the dominant mode of energy loss of a cosmic string network is particle 
production and not gravitational radiation as the loops collapse almost immediately. 
Indeed, it may not be possible to talk about loop production at all: the loop production 
we observe favours loops formed at cusps, which could annihilate into particles before 
the loop is formed |15| . Furthermore, gravitational radiation directly from a network 
without a small scale ( -C £ is negligible. 

Recent calculations |17| on the possible role of GUT quanta decaying into high energy 
cosmic rays assumes that particle production is sufedominant to gravitational radiation 
as a means of energy loss, and conclude that GUT strings could not be an appreciable 
source of cosmic rays. In our suggested energy loss scenario, we may come up against an 
observational bound of high energy cosmic rays. Simple models of GUT quanta decay 
processes give the upper bound for the proportion of the string energy going into high 
energy cosmic rays as O(10 -3 ) or less [lfj. Such calculations may put a significant 
constraint on GUT models if the GUT scale string scenario is to be viable. 
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